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Abstract. The properties of dynamical Dark Energy (DE) and, in particular, the 
possibility that it can form or contribute to stable inhomogeneities, have been widely 
debated in recent literature, also in association to a possible coupling between DE and 
Dark Matter (DM). In order to clarify this issue, in this paper we present a general 
framework for the study of the nonlinear phases of structure formation, showing the 
equivalence between two possible descriptions of DE: a scalar field cj) self-interacting 
through a potential V((f>) and a perfect fluid with an assigned negative equation of 
state w(a). This enables us to show that, in the presence of coupling, the mass of DE 
quanta may increase where large DM condensations are present, so that also DE may 
partake to the clustering process. 
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1. Introduction 

An important part of cosmology concerns the large-scale clustering of matter in galaxies 
and galaxy clusters and its statistics. The study of Large Scale Structure (LSS) can 
reveal many aspects of the physics of the early Universe as well as its matter content 
during the cosmic history. 

Several observations made over recent years, related to a large extension to LSS and 
anisotropies of the Cosmic Microwave Background (CMB) as well as the magnitude- 
redshift relation for type la Supernovae pQ, have given us a convincing picture of the 
energy and matter density in the Universe. 

Baryonic matter accounts for no more than 30% of the mass in galaxy clusters while 
the existence of a large clustered component of Dark Matter (DM) seems now firmly 
established, although its nature is still unknown. However, they contribute to the total 
energy density of the Universe with only a few percent and about 25% respectively. 

No more than another few percent could be accounted for by massive neutrinos, 
but only in the most favorable, but unlikely case. According to [2] the total mass of 
neutrinos cannot exceed the limit of 1.43 eV. A very small part (1CT 4 ) of the total 
energy density is due to massless neutrinos and CMB radiation. 

The model suggested by observations is only viable if the remaining 75% is ascribed 
to the so-called Dark Energy (DE) responsible for the present day cosmic acceleration. 

Although strongly indicated by the observations, the existence of DE is even 
more puzzling than DM. It can be identified with a cosmological constant or with a 
yet unknown dynamical component with negative pressure. On the other hand, its 
manifestation can be interpreted as a geometrical property of the gravity on large scales 
denoting a failure of General Relativity (GR) on those scales (see [3] for a review). 

Within the context of GR, as alternative to the cosmological constant, DE is usually 
described as a self-interacting scalar field or a cosmic fluid with negative pressure (see 
[I] and references therein). It is usually assumed that density perturbations of DE play 
a negligible role in the structure formation because of its very small mass ~ H (H being 
the Hubble parameter). Accordingly, perturbations should appear only on very large 
scales and are bound to be linear so that rates of structure formation and their growth 
are influenced by DE only through the overall cosmic expansion [5]. 

Nevertheless, the clustering properties of dynamical DE have been subject to recent 
debate. The behaviour of DE in the presence of high concentration of matter and in 
particular of DE-matter coupling [HI [7] , is not clear. 

Then, the key question is whether DE actively participates in the clustering and 
virialization processes developing non-linearity on relevant scales. Some attempts to 
solve the problem have been done in [HI [9] and an analytical expression for uncoupled 
DE density perturbations valid both in the linear and nonlinear regime was derived in 

ma. 

The aim of this paper is to provide a general framework which can be useful for the 
study of the nonlinear phases of structure formation in cosmologies where a coupling 
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between DE and matter is present. The paper is organized as follows. 

We assume a Friedmann Robertson Walker (FRW) background Universe. Because 
we are interested in the epoch when structures form, we neglect the radiation and 
consider a system of two coupled fluids. A pressureless fluid and one with pressure 
which can be interpreted as DM and DE respectively. Equations for the background 
evolution are provided in Sec. [2j Then, in Sec. [3X starting from a general relativistic 
treatment, we derive the Newtonian limit of the nonlinear perturbation equations for the 
density contrast. In sec. HI following [TT] (see also [12]), we review the scalar field-fluid 
correspondence showing that the formalism described in the previous sections apply also 
to a scalar field. It is then applied to the coupled DE model in Sec. El We conclude in 
Sec. El 

2. Background equations 

Assuming GR holds, let us start from the background equations. We adopt the following 
conventions: the four spacetime coordinates are labelled by Greek letters running from 
to 3 while Latin letters run from 1 to 3, labelling the spatial coordinates. We are 
interested in the epoch when structures form so that radiation is neglected. We consider 
a Universe filled with baryons (b) and two coupled perfect fluids with vanishing and non- 
vanishing pressure which, in case, can be interpreted as DM and DE respectively and 
are labelled with dm and de- According to general covariance, the sum of the individual 
stress energy tensors T^ v (a = b, dm, de) must be conserved so we can write: 

^ t^{b)u = ^ \^{dm)v = ~Cv V \T[de)v = C v (1) 

Here V M is the covariant derivative while the 4-vector C M parametrizes a possible 
interaction between DM and DE being the energy and momentum exchanges described 
by Co and C, respectively. For each component, the stress-energy tensor can be put 
in the perfect fluid form: 

TZ = p5Z + (p + p)u»u u (2) 

Here p and p are the energy density and pressure of the fluid as measured in the 
rest frame and w M is its 4-velocity. We also assume the coupling term C v = 
C v ({vi, v n }dm, {vi, v m }de) to depend on the dynamical variables of the two coupled 
fluids that we have generically indicated with v a (a = 1, 2, ...). 

In a spatially flat FRW background with metric ds 2 = a 2 (r) (—dr 2 + dx l dxi) (r is 
the conformal time) the Friedmann equation for the scale factor a and equations (JT]) 
read: 

ft 2 = ^-G (p b + p dm + Pde) a 2 (3) 

p b + 3Hp b = p dm + 37-lpdm = -C pde + 3ft [fide + Pde) = C (4) 

where ' denotes the derivative with respect to r, H = a/a and ~ denotes the background 
quantities. For a generic function F(v i, v 2 , ■■■) then we have F = F(vi,v 2 , ...) so that 
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C = C ({v\, ...}dm, {vi, ■■■}de)- Notice that from (JTJ) and fl2]) it follows that the coupling 
term in the FRW background is C* M = (Co, 0, 0, 0). 

3. Perturbation equations 

According to the standard picture, the observed structures have formed via gravitational 
instability from the tiny inhomogeneities left by cosmic inflation. 

Relativistic theory of cosmological perturbations [T3], |H] provides a useful tool to 
understand the evolution of inhomogeneities relating the physics of the inflation to CMB 
anisotropy and LSS. 

Despite relativistic effects, the physics of the evolution of perturbations is quite 
simple during the period from the end of inflation to the beginning of nonlinear 
gravitational collapse as they have small amplitudes and different wavelengths evolve 
independently. 

Anyway, linear theory breaks down when matter inhomogeneities become 
significantly denser than the cosmic medium. Although the Newtonian theory is 
appropriate to describe the dynamics on scales small compared to the Hubble radius 
(but larger than the Schwarzschild radius), full numerical simulations are needed to 
follow the nonlinear phases of structure formation (a comprehensive review on numerical 
techniques is provided, e.g, by [15]). Only in cases where some symmetries are present, 
analytic or semi-analytic schemes apply. The simplest approach is to follow a spherical 
density enhancement; although real fluctuations are expected to be highly irregular and 
random, much work has been done along this line, starting with [161 Ej who studied 
the problem within the frame of pure cold DM models. Then, the results have been 
generalized to the case of ACDM [18J and other DE models [5J El fl9] . In spite of the 
lack of cogent physical motivations, spherical geometry yields results suitable to work 
out useful mass functions and to study their evolutions within the Press & Schechter 
[20] - or similar [21 J - approach (see, e.g, [22]). 

The aim of this section is to provide a general framework which can be used to 
study the nonlinear phases of structure formation in the context of coupled DM-DE 
models. Without considering any specific coupling between the dark components we 
derive and discuss a set of equations that help to understand the role of the DE and its 
possible coupling with DM in the clustering process. 

Dealing with LSS it is sufficient to consider weak gravitational fields and 
nonrelativistic gravitational sources. Note that, because of the weak fields, slow motions 
do not necessarily imply small density fluctuations. 

Starting from a general relativistic treatment, we derive now the perturbation 
equations in the Newtonian limit, valid to all orders in the density contrast. We choose to 
work in the conformal Newtonian gauge since it most closely corresponds to Newtonian 
gravity and consider only scalar perturbations. In this gauge, the metric of a perturbed 
flat FRW Universe reads: 



ds 



2 




2 



2 



(5) 
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where $ plays the role of the Newtonian potential, \I/ is the Newtonian spatial curvature 
and |$|, |^| << 1. 

Consider now a generic perfect fluid. In the nonrelativistic limit its coordinate 
velocity v l = dx l /dr can be treated as a perturbation of the same order as the metric 
perturbations. Then, keeping only terms linear in and \1/ and neglecting terms 
involving products of v % and the metric potentials the components of the 4-velocity 
= dx 11 / \/—ds 2 become: 

u = a-\l-$) u i = a~ 1 v i u = -a(l + $) u i = av i (6) 

With the same approximations, we can now write the components of the stress-energy 
tensor (J2j): 

T ° = -p T?=(p + P y Ti = -(p + py T]=p§) (7) 

Here p = p(r) + Sp(r, x l ) and p = p(r) + <5p(r, x % ) are the perturbed energy density and 
pressure of the fluid. Notice that we have not required 5p and 5p to be small compared 
to the background quantities. 

The equations of motion follow from [17] : 

= I (d u g^)T^ - -L d, L/=g Tfj = C u (8) 



2 V^g 

where g is the determinant of the metric g^ v given by ([5]) and to the first order in 
the gravitational potentials y/—g = a 4 (l — 3\l/ + $). Possible interactions with other 
components are described by C M = C M (r)+5C Ai (r, x l ) where 5C represents deviations from 
the value C taken in an homogeneous and isotropic Universe . Taking C M = 0, C M , — C M 
equations (OQ) are then recovered. From the component v = of (jHJ) it follows: 

% + 3 (h - ™) ( P + p) + V • [{p + p)v] = C (9) 

while the component v = % gives: 
d 

tt- (p + p) v ] + w (p + p) v + Vp + (p + p) V$ = -C (10) 
or 

( | v | 2 = fjt> 1 , |C| 2 = CiC 1 ). Then, making use of the unperturbed part of (jHJ): 

p-+3-(p + p)=C (11) 
a 

together with the relation p + p = p(l + uj) + 5p(l + cj p ), where lj = p/p is the equation 
of state parameter and co> p = Sp/5p, we can extract from ([9]) and ( TTUT) the equations for 
the density contrast 5 = Sp/p and the peculiar velocity v: 

6 + 3H(u p - u)8 - 3*(1 + R5) + (1 + cu)V • [(1 + R5)v] = C °~ C °_ (1 + ^ 

L.. . 1 f(Cn -Cn)R + C n UJ+uJ _xM 

V + 



H(1 - 3M) + _i_ (Co-C»)fi + C + ^ + ^ 
' 1 + R5 \ p 1 + u . 



v 



+ Y*R + v$ = % (12) 
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Here R = (1 + ui p )/{l + oS) and V denotes the spatial gradient. 

Notice that in the small perturbation limit (dp/p, 5p/p << 1), the parameter uj p 
is related to the sound speed c 2 s in the fluid. In general, Sp = Sp a + Sp na where the two 
terms on the r.h.s are the adiabatic and non-adiabatic contributions to the pressure 
perturbations and c 2 s is then their propagation speed in the fluid rest frame (rf), i.e. 
the frame where TP = 0: 



r 2_ *E 



(13) 

rf 

After performing a gauge transformation from the rest frame gauge to the conformal 
Newtonian gauge one arrives to the following relation between 5p and 5p [T4"t [231 El] 
(we assume c 2 s to depend only by r): 

V6p = V6 Pa + V6 Pna = c 2 s V6p + (c 2 s - c 2 J[V8p + pv] (14) 
which can be rewritten as: 

V5p = c 2 s V5p - (c 2 - c 2 J[3H(l + u)p - C ]v (15) 
showing that the coupling enters Sp explicitly. In the above relation, 

c 2 = t = uj + 6ji (16) 
° P P 

is the adiabatic sound speed which coincides with c 2 if the fluid is barotropic (p = u(p)p). 
Pure adiabatic perturbations, however, go unstable for fluids with uj < 0, e.g. DE 
fluids. Therefore, in order to ensure their stability (c 2 > 0) and have a DE component 
phenomenologically viable, non-adiabatic modes must be also present (notice that the 
condition c 2 > needs to be imposed by hand if DE is assumed to be a fluid while 
Cg — 1 follows without assumptions in the case of DE scalar fields. Also notice that 
a scalar field has in general a non-barotropic equation of state (see [241 [25] and next 
section)). Nonetheless, non-adiabatic instabilities could arise in the presence of strong 
DM-DE couplings or on super-Hubble scales. Adiabatic and non-adiabatic instabilities 
have been discussed in [26J. 

Equations (|12[) can be generalized to the case in which anisotropic stresses are 
present, i.e imperfect fluids. In that case we have T ? = pSj + II*- where II* is the 
anisotropic stress tensor which satisfies 11^ = and U^u u = 0. The equation for the 
peculiar velocity evolution must then be modified adding the term VII to the l.h.s. 
(see [1 3j for further details). In the next section we show that a scalar field minimally 
coupled to gravity has vanishing anisotropic stresses and its stress energy tensor can be 
put in the perfect fluid form. It follows that a DE scalar field can be described as a 
perfect fluid. On the other hand, if one assumes DE to be a fluid, anisotropic stresses 
can not be a priori excluded. Anyway, as we shall see shortly, the absence of anisotropic 
stresses leads $ = ^. 

In additions to (fl2"]) . one has to consider the equations of metric perturbations 
obtained from the perturbed part of the Einstein's equations. Because of the weak 
field approximation, the Einstein's tensor will be only linearly perturbed while 
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the perturbations in the stress energy tensor components of each fluid a follow by 
subtracting their background values Tjf from those obtained from (jTJ): 

SG^SnG^iT^-f^) (17) 

a 

Einstein's equations then read: 



V 2 (* + W$) = -47tGa 2 di 



y + H(<3? + 2tf ) + (2- - H 2 ] $ - -V 2 (^ - $) = — Ga 2 V <fy Q 
V a J 3 3 ^ 

<9^(* - $) = 87rGa 2 £ rf (a)j (18) 

Thus, if anisotropic stresses are absent the last equation shows that \1/ = $. In this 
case, focusing on scales well below the Hubble radius, we can neglect the second term 
on the l.h.s. in the first equation, because 7^$ ~ 7^ 2 $ << V 2 $. This correspond to the 
Newtonian limit and the first of equations (fTBl reduces to the usual Poisson equation: 

V 2 $ = 47rG]T<fy Q (19) 

a 

In our particular case, neglecting the term proportional to \l/ in the first of ([121) 
because of the Newtonian approximation and replacing the term 5 in the second, the 
perturbation equations for our three perfect fluid components are: 



Baryons (lu m 0, c s negligible after recombination, = 0) 

h + v ■ [(1 + s b )v b ] = 

v b + Hv b + V$ = (20) 

DM (Pressureless perfect fluid) (to = 0, tu p = 0, C M = — C M ) 

C -C (l+5 

dm) 



$dm + V • [(1 + 5 dm )v dn 



Pdn 



+ ( n ~ - rrh s ) Vdm + v$ = - rrr, \ ( 21 ) 

V Pdm(l + Odm) ) Pdm{i-+0dm) 
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DE (Perfect fluid with pressure) (u = p de /p de , u p = Sp de /5p de , = C M ) 



5 de + 3H(u p - u)5 de + (1 + oj) V ■ [(1 + R6 de )v de ] 



Cq — Co(l + &dm) 



Pde 



1 — 3u + (1 — 3uj p )R5 de 
1 + Mfe 



+ 



(l+u)(l + R6 de ) 




+ 



(Cq — Cq)-R + C 

p de (l + RS de ) 



+ V$ 



c 



(22) 



p de (l+w)(l + J R5 de ) 



p de (l + uj)(l + R5 de ) 



Equations (|20p . (|2ip and (|22p are valid to nonlinear order in the density contrast since 
we have not assumed 5 << 1. They can be used to follow the evolution of a collapsing 
region in nonlinear regime and could be useful to study the clustering properties of DE 
by using the spherical collapse approach. The problem of spherical fluctuation growth 
in coupled DE models has been considered in [19] . There, equations ( 1201) and ( 12TT) were 
used to follow the evolution of a set of concentric shells of a spherical top-hat overdensity 
in baryons and DM while the DE was kept homogeneous. The shell description of the 
overdensity was needed because the coupling causes DM particles to have a different 
dynamics than the baryons. This can be seen from an inspection of the above equations. 
The most significant result of that analysis is the high level of segregation between DM 
and baryons arising during the growth of fluctuations: up to ~ 60% of baryons can be 
expelled from the fluctuation in halo encounters before its final virialization. 

The same approach could be used relaxing the assumption of homogeneity for DE 
and including in the system an additional set of concentric shells for DE. Their evolution 
is then described by ( 1221) . In contrast to DM, dynamics of DE shells is affected, not only 
by the coupling but also by the non-vanishing pressure of DE itself. We shall consider 
this problem in a subsequent paper. 

A similar set of equations was used in [27] to study nonlinearities in the case of 
uncoupled DE . However, the authors consider the pseudo-Newtonian approach of [28J 
and the DE equations they use lack of some terms with respect to ( 1221) . 

In the next section we show that that equations (|22|) applies equally well when DE 
is described by a scalar field. 

4. Scalar field 

In most DE models, DE is due to a scalar field either minimally or non-minimally 
coupled to matter and/or gravity. In this paper we do not discuss coupling to gravity, 
while coupling with matter will be considered only in a subsequent section. Anyway, 
results of this section are not affected by the presence of this latter kind of coupling and 
are fully general. 

Here we show that the equations deduced in the previous section are valid also for 
a scalar field. In particular, following [IT] , we show that the stress energy tensor of a 
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scalar field, if minimally coupled to gravity, can be put in the perfect fluid form. As a 
consequence equations (|T2|) as well as (1221) can be applied. 

Let us start with a minimally coupled scalar field self-interacting through a 
potential V((f>). In a spacetime described by a metric g^ v with signature ( — h ++) its 
Lagrangian reads: 



C 



-g 



-V(<f>) 



and the stress energy tensor takes the form: 



T? = d^(j>d v (j> - 5» 



+ V 



(23) 



(24) 



The equations of motion then follow either from V M X^ 
equation: 

1 



-A, 



-g 



-gd" 



- V = 



or from the Euler-Lagrange 



(25) 



(the prime denotes the derivative with respect to 0). 

Assuming that the vector <9 M is timelike, namely d^ipd^cj) < 0, one can assign a 
4-velocity u^ L to the scalar field. It can be chosen as the unique timelike vector with 
unit magnitude constructed from d^cf) (for more details see [Tl~[ [12]): 

d»4> 



1 



(26) 



From this choice of the 4-velocity, it follows that the scalar field <fi can be represented 
as a perfect fluid and the stress energy tensor ([24]) can be put in the form ([2]). In order 
to show this, let us start from the stress energy tensor of a general fluid: 



Tp» = pg^v + (p + p)u^u u + q^Uv + u^q v + £ 



[IV 



(27) 



where possible anisotropic stresses are accounted for by the heat flux vector q^ and the 

i that for a 



tensor £^ which satisfy = 0, = and T,^u u = 0. We show now that for a 



scalar field these terms vanish. Making use of the projection tensor 
one can extract from ([27]) all the relevant quantities: 

1 



p = T ilv u iL u v 



V 



-T uu h^h u a 



T aX h;h x u - P h 



By using ([241) and ([26]) . the above quantities can be found for the scalar field 4>: 

1 



P 



+ V 



V 



2 9 » 



V 











(2f 



(29) 



Both the heat flux and the anisotropic stress tensor vanish, showing that the scalar field 
stress energy tensor has the usual perfect fluid form. 

However, it is important to mention that, although this mathematical equivalence 
relates a scalar field with perfect fluid models, differences remain in their physical 
interpretation. In general, a fluid is a continuum model deduced from thermodynamical 
considerations of a system containing a significant number of particles interacting 
through elastic collisions. Its continuum variables (e.g., p, u^) are defined in terms of 
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fundamental discrete variables in a limit procedure. Once the large particle number limit 
is violated the model is expected to fail. Because of the thermodynamical motivation, 
physical constraints on p and w M are then imposed, namely that p > and that 
is timelike. On the other hand, a scalar field model is a continuum model and the 
solutions of its equations are generally not subjected to the constraints above mentioned. 
Thus, the assumption d^cpd^cp < only refers to those scalar field solutions that can be 
interpreted as physical fluids. 

Let us now consider a perturbed FRW Universe described by the metric (JSJ). 
Equation of motion (|23|) then reads: 

4> + (2H - $ - 3^) - (1 + 2^ + 2$)V 2 + a 2 (l + 2$)V = V($ - ■ V0 (30) 
and the components of the fully nonlinear stress energy tensor are: 

A2 



± 



2a2 (l-2*) + ^(l + 2*)-V 



2? = -4(1 -2$)$. 



1 + 2tf ) + <J* 



^ (1 - M) -5Af (1+2t) -v 



+ 2-0 + aV = 
a 



(31) 



2a? y ' 2a? 

Then, decomposing <fi as the sum of an unperturbed part and a perturbed one 5$: 

= 0(r)+50(r,x i ) (32) 
the Oth-order terms of (1301) give the equation for the homogeneous part of 4>: 



(33) 



while, assuming no anisotropic stresses ($ = \1/), the higher order terms yield for the 
perturbations: 

5(j) + (2H - 4$) 8(j) - 4<±>0 - (1 + 4$)V 2 50 + a 2 \(l + 2$)V* - V"' 







(34) 



From (131]) one can also work out the background energy density and pressure and the 
corresponding perturbations: 

$ 2 - $ 2 is'd), v 5(f) 2 , x J2i(di5d)) 2 , 



2a? 



2a 2 



^^~ v * = -^ + ^F(i-») + ^(i-»)- 5 ^(i + »)-n' 

(35) 

where 5V = V — V" and V" = V(0). The scalar field peculiar velocity v % is obtained by 
comparing the expression (1261) . which can be rewritten as: 



vJ 1 = a' 1 



;i + 2*)-(i-2*)s(|Q" 



(36) 
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Table 1. Scalar field-perfect fluid correspondence in the small velocity and weak 
gravitational field limits 

Perfect fluid Scalar field 



'<t>„-i 



dr 



t)",p 



-P 

-(p + p)v l 



-^(1-2$)-^) 



^(1-2$)-^) 



with the fluid 4-velocity: 

dx^ 



vr 



-i 



Taking 



z, it then follows 

{ dx % d\ 
dr 



2$) 




dx^ 
~dr 



d°(j) 



(37) 



(3? 



As we are interested in the small scale and weak gravitational field limits some 
approximations are possible. Terms proportional to $ can be dropped in equation (134"|) . 
We can also neglect the time derivatives of the gravitational potential with respect to 
7i, since $ ~ $7i. We are then left with: 

5(f) + 2H5(f) - V 2 5(f) + a 2 5V = 4$0 (39) 

The correspondence between the perfect fluid and scalar field descriptions is 
summarized in the table [1] where the expressions taken by the 4-velocity and the stress 
energy tensor components in the small velocity and weak gravitational field limits are 
given. Notice that these limits, in the scalar field description, correspond to considering 
only linear terms in <9 M 0/<9°0, \1/ and $ and neglecting terms involving products of the 
former and the metric perturbations. It is worth to mentioning that equation (|39l) is 



equivalent to the first of ( 1221) . from which it can be directly derived by setting C M = 0. 
However, once the prescriptions of Tabled] are taken into account, equation (1391) is more 
simply deduced from (jSJ) or directly from after considering the small scale limit and 
eliminating the homogeneous part. 

As mentioned in the previous section, in general, a scalar field cannot be considered 
as a barotropic fluid being its equation of state given by: 

p = p - 2V (40) 
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Therefore, for small perturbations, if V ^ const, we have 5p/5p ^ p/p = c 2 yielding 
also non-adiabatic contributions to the fluctuations. The sound speed is given again 
by ( |T3i) where the rest frame is defined by the hypersurfaces = const orthogonal to 
the rest frame 4- velocity (1261) . By definition, in the rest frame, the scalar field carries 
no perturbations (S(p\ r j = 0) leading SV\ r j = 0. It follows that energy density and 
pressure perturbations come purely from kinetic terms so that o~p\ r f = &p\ r t an d c 2 s — 1 
independent of the form of V(<p) (for a detailed discussion see [29]). 

In the next section we introduce the coupling through a discussion on the coupled 
DE model. 

5. Coupled Dark Energy 

In addition to self-interaction, a scalar field can in principle be coupled to any other 
field present in nature. However, if this field is the one which accounts for DE, in order 
to drive the cosmic acceleration, its present time mass is expected to be, at least on 
large scales, ~ Hq ~ 10 _33 eK. Such a tiny mass gives rise to long-range interactions 
which could be tested with fifth-force type experiments. Couplings to ordinary particles 
are strongly constrained by such a kind of experiments but limits on the DM coupling are 
looser (constraints on coupling for specific models were obtained in [30i [3TJ from CMB, 
N-body simulations and matter power spectrum analysis). If present, DM coupling could 
have a relevant role in the cosmological evolution affecting not only the overall cosmic 
expansion but also modifying the DM particles dynamics with relevant consequences on 
the growth of the density perturbations in both linear and nonlinear regime (e.g., on 
halo density profiles, mass function and its evolution) [3QI [32j [33] . Here we review one 
of the most popular models where a coupling between DM and DE is present, namely 
coupled DE [32], and show how the formalism introduced in Sec. [3] can be applied. 

The model is quite general and includes a vast class of models also motivated by 
string theory as proposed in [6l [7] . Conservation equations for baryons, DM and a DE 
scalar field <j) read: 



= V„7f dw> = -c(0)T (dm) V,0 = c(0)T (dm) V,0 (41) 



where T^ m ) is the trace of the DM stress energy tensor and c(0) parametrizes the 
strength of the DM-DE interaction. Such a coupling can be derived from a scalar- 
tensor theory after performing a conformal transformation (see [3l] for a review on 
scalar-tensor gravity). Different couplings were also considered by various authors [35J. 

The equation of motion of the scalar field follows from the last of (j4"Tj) and fl24j) or 
equivalently from the lagrangian 



so that the coupling C M introduced in Sec. [2] now takes the form: 
c » = c (0)?W)V M 



(42) 




(43) 
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in the case when c(0) = d(lnB) / d<f>. The last term in (14"3"|) sets the coupling between the 
DE and the spinor field \ supposed here to yield DM (results are however independent of 
whether the DM particles are scalars or fermions). Notice that in such a theory the mass 
of the DM particles, m dm = £>(0)m, as well as their energy density pd m = m dm < XX > 
are 0-dependent. In this case equation ( 1391) reads: 

5xj) + 2H5(j) - V 2 5(f) + a 2 b~V' eSS = 4$0 (44) 

where 5V e ff = V e ff — V e ff and V e ff = V + pd m . Notice that no assumptions obout the 
value of 5cj) were made in deriving the above equation which is valid also for nonlinear 
5<j). 

It is also worth to mentioning that, unlike the uncoupled case, the effective mass of 
the scalar field, = V"^{(j)) = V"(<p) + d(<f))pdm + c 2 ((j))pdm, as well as the value of </> 
depend on the local DM density pd m . It follows that, in massive objects, where the DM 
density is high compared to the background density, could be very different from 
the value it takes in the cosmos where pd m ~ p cr (p cr being the critical energy density) 
and <$ = <$. Thus, perturbations in the DE energy density and pressure could be non- 
negligible on the scales relevant for the structure formation. Effects of a scale-dependent 
mass was considered in the contest of the so-called chameleon theory [36] where, unlike 
the coupled DE case, the scalar field universally couples to all the kinds of matter, e.g. 
DM and baryons. However, even though this theory simultaneously provides a viable 
cosmology and a mechanism which allows the scalar field to evade constraints from fifth 
force effects, a very low value of the energy scale entering the chameleon potential is 
needed. Although tracker potentials were considered, it seems therefore hard to avoid 
fine-tuning and coincidence problems. 

Alternatively to the equation (jHJ), in order to study the evolution of DE 
perturbations, one can use the perfect fluid description and the equations ( 1221) once 
the coupling vector C M is rewritten in term of p de , Pde and v de . Let us consider first the 
case when the parameter c of ( 1421) is constant. In this case, by using the relations of 
table [U the coupling vector becomes: 

C = Cp dm (l + 5 d m)[pde(l + W)(l + Rb~de)]*CL 



C = -Cp dm (l + 5dm)[pde(l + W)(l + RSdej^dVfe (45) 

Problems can arise when c is dependent. However, if the potential V(4>) is a 
strictly monotonic function, as in the case of an inverse power law or an exponential 
potential, one can invert it to obtain <fi as function of V = (pd e —Pde)/2- As an example, 
consider the Ratra-Peebles potential V = A 4+a /<j) a where a > 0. In this case, 

/A 4+a \" _ / 2A 4+Q 

~ \pde[{l-u) + 5 de {l-u p )} 

so that c = c(pde — Pde] A, a). Notice that, once the energy scale A is fixed, the a 
parameter follows by requiring p rfe (or p m = pdm,o + Pb,o) to have a value as inferred 
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by observations (here the subscript o refers to the present time). Alternatively, one can 
fix a and f>de,o (or p m; o) then find A. 

6. Conclusions 

One of the main topics which has attracted a deep interest in recent years concerns 
the clustering properties of dynamical DE, in particular when a coupling with matter 
is present. Scalar fields are among the favorite candidates for DE. However, in certain 
circumstances, it can be more advantageous describing DE in terms of a perfect fluid 
rather than in terms of scalar fields. 

This paper is a technical review meant to set up a formalism for the study of 
nonlinearity in the presence of DM-DE coupling. To this aim, we derive the Newtonian 
limit of the perturbation equations for DM and DE (and baryons) valid to nonlinear 
order in the density contrast, highlighting the equivalence between the perfect fluid and 
the scalar field descriptions of DE. We then provide some prescriptions (summarized in 
table d]) to relate the two descriptions and apply them to the coupled DE case. This 
enables us to give the expression of the coupling vector in terms of the fluid variables 
and show that, in the presence of coupling, the mass of DE quanta could increase where 
large DM condensations are present. As a consequence, also DE can take part in the 
clustering process. 

The framework here presented gives the DE equations of motion in the presence 
of coupling once the time evolution of the DE equation of state, 00(a), is assigned, 
without reconstructing the self-interacting potential V(4>). It can be useful to follow 
the evolution of a collapsing region in nonlinear regime in order to study the clustering 
properties of DE, e.g. by using the spherical collapse approach. This point is currently 
under investigation. 
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